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Abstract 

We discuss general properties of D-brane model building at toric singularities. 
Using dimer techniques to obtain the gauge theory from the structure of the singu- 
larity, we extract results on the matter sector and superpotential of the correspond- 
ing gauge theory. We show that the number of families in toric phases is always 
less than or equal to three, with a unique exception being the zeroth Hirzebruch 
surface. With the physical input of three generations we find that the lightest 
family of quarks is massless and the masses of the other two can be hierarchically 
separated. We compute the CKM matrix for explicit models in this setting and find 
the singularities possess sufficient structure to allow for realistic mixing between 
generations and CP violation. 
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1 Introduction 



Explaining the structure of the Standard Model of particle physics is one of the biggest 
outstanding problems in string phenomenology. The success story of the Standard Model 
is uncontested but an ultraviolet completion, as in principle provided by string theory, 
is an important challenge. 

There are many different proposals for embedding the Standard Model and its exten- 
sions into various classes of string models. Besides the general features of the Standard 
Model, such as the gauge symmetries and spectrum, there are also many observables 
that should be calculable from the corresponding string model and confronted with ex- 
perimental results. In particular, realistic models must explain the hierarchy of masses 
of the different generations of quarks and leptons and the information encoded in the 
CKM matrix. 

Efforts in string model building go back to heterotic string theory and include more 
recently constructions in brane engineering within type II flux compactifications. An 
attractive feature of type IIB model building is the presence of successful scenarios for 
moduli stabilisation [TJ[2] . The LARGE volume scenario (LVS) [2] is particularly attrac- 
tive: the LARGE volume ensures that perturbative corrections are well under control 
and does not require fine-tuning of the flux superpotential. The hierarchy problem can 
be solved by gravity/moduli mediated supersymmetry breaking [3--6J. Depending on the 
overall volume of the bulk geometry one can lower the string scale to values smaller 
than or similar to the GUT scale Mqut ~ 10 16 GeV, making them attractive from the 
perspective of gauge unification. 

The existence of a large volume implies that the Standard Model has to be located 
at a local region within the full Calabi-Yau space [TJ. It is natural to think of those 
models as supersymmetric models, with supersymmetry breaking dealt with from the 
bulk [6,8,9J. We hence will, from now on, talk about supersymmetric model building 
and expect supersymmetry to be broken by bulk/moduli effects. 

This viewpoint is applicable to all model building efforts in this category: intersecting 
brane models [T0] - [T3] . F-theory models [TH422] and branes at singularities [7tl8| [23l427] . 
Several models with realistic matter content and potentially couplings have been devel- 
oped in these scenarios. However, obtaining the correct flavour physics remains challeng- 
ing to date. For example in the widely studied F-theory models, the flavour structure is 
not complete at tree-level [28H33] and one has to introduce non-commutative fluxes [3l] 
or appropriate instanton effects [35] to achieve for Yukawa matrices with rank greater 
than one. 

In this article we concentrate on the class of models corresponding to D-branes at 
singularities which represents one of the most promising avenues to obtain the Standard 
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Model from string theory: 

• The models are truly local in the sense that many features can be addressed in a 
way that is essentially independent of the moduli stabilisation problem. 

• Even though the distances involved are in a regime smaller than the string scale, 
string theory is well under control. 

• In the case of IIB string theory, the models can be determined from local F- 
theory models in the limit of vanishing cycles, yet their properties are much more 
constrained and in some sense represent a minimal class of F-theory models. In 
particular there is much less freedom to include the matter sector and couplings, 
making realistic model-building more rigid. If a realistic property is obtained, its 
presence can be argued to be more robust than in local F-theory models where 
there is much freedom to engineer the models. For instance simple GUT models 
are not possible to obtain with non- vanishing quark masses. However, several other 
generalisations of the Standard Model are easy to obtain. Contrary to simple group 
GUT models, Yukawa couplings giving rise to fermion masses are not forbidden 
by anomalous £7(1) symmetries and the structure of these couplings is elegantly 
determined in terms of the corresponding dimer diagrams [36til2] . 

• In principle, if the standard model is on D3 branes at a singularity, gauge coupling 
unification is automatic since, to leading order, all the gauge couplings of a product 
gauge group are given by the expectation value of the dilaton field, a property 
shared by heterotic models but not by intersecting brane models. The question 
of unification reduces to whether the matter sector is such that the measured low 
energy values of the gauge couplings can be obtained after renormalisation group 
running. Notice that a simple GUT group is not needed to achieve unification. 

• To leading order a vanishing value of the blow-up mode is a natural extremum of 
the effective action and then considering the model at a singularity is a promising 
ansatz once moduli stabilisation is taken into account. 

• Over the past few years very powerful techniques have been developed [36-42J to 
describe these models in terms of quiver and dimer diagrams that provide the 
relevant information concerning the spectrum and couplings of the corresponding 
gauge theory. These techniques to study branes at general singular points have been 
used to address several aspects of the AdS/CFT correspondence as well as local 
supersymmetry breaking and properties of M2-branes, but, except for the simplest 
cases of orbifold singularities and some del Pezzo surfaces, their phenomenological 
aspects for string model building have not been explored. 
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• The gauge theories that can be obtained at the singularities are highly restricted. 
Both the matter content and superpotential are completely determined. One may 
only vary the ranks of the gauge groups and vevs of fields (corresponding to reso- 
lution of cycles). 

1.1 Summary of Results 

In this paper we focus on local model building with branes at singularities and de- 
scribe general constraints on flavour physics within a large class of these models, the 
so-called toric singularities for which the matter content and superpotential can be ob- 
tained systematically. We start with a detailed introduction to the relevant aspects of 
toric singularities in terms of quiver and dimer diagrams. In particular we describe how 
to determine the matter spectrum and perturbative superpotential for D3 branes at the 
singularity, including also D3-D7 states. The most efficient way to determine the gauge 
theory from the toric singularity is through an algorithm [43j which we describe in detail 
since it shall be the primary tool used in obtaining our results. Our main results can be 
summarised as follows: 

Three families bound. We find as a general result that the number of families in 
toric models is bounded by Nf = 3. The only exception is the well known case 
of the zeroth Hirzebruch surface which has four families, but has a Seiberg dual 
phase with two families^] The bound Nf < 3 was found before in singularities, 
with Z 3 saturating the bound [7]. It is remarkable that this result extends to the 
much larger class of toric singularities. Given the arbitrariness in the number of 
families in most string constructions, it is intriguing that the physical value plays 
an important role in this class of models. In principle, non-toric phases obtained 
by Higgsing could generate models with more families. 

Hierarchy of masses. The explicit knowledge of the superpotential allows us to 
compute the mass matrix for the quarks. When a physically realistic choice of 
quarks is made, we find that there is always one vanishing mass eigenvalue in a 
toric singularity. Generically it is possible to find hierarchical masses for the other 
two eigenstates (except for the zeroth del Pezzo singularity dP ). For three families 
we find then the masses (M, m, 0) with M » m. This result was found for dP\ 
in [8j and it is remarkable that extends to the general toric case. 

CKM matrix. We have also computed the CKM matrix for two classes of models at 
toric del Pezzo singularities: the first is where both the up and down quarks arise 

: As the example of the zeroth Hirzebruch surface suggests the origin of the 3-families is not due to 
the 3-complex dimensions of the bulk geometry. 
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from D3-D3 states, and the second has one type of quarks as D3-D3 states and the 
other type as D3-D7 states. We can construct the correct CKM matrix in both 
type of models by appropriate values in the ratios of the Higgs vevs. We illustrate 
our findings with concrete examples based on SM-like and left-right symmetric 
models at del Pezzo singularities. We show that in the first class of models the 
dPi singularity allows for the correct flavour mixing. In the second class of models 
we find the correct mixings for the (IP2 and dP% singularities. In our analysis we 
neglect perturbative corrections to kinetic terms or non-perturbative corrections 
to the superpotential which in principle could influence the small mixing angles in 
the CKM-matrix. 

CP violation. Given the structure of the Yukawa matrices, the amount of CP violation 
can be easily computed in terms of the Jarlskog invariant J [44,45j. With the 
hierarchical structure in the CKM matrix, the magnitude of the Jarlskog invariant 
is automatically in the desired range [16]. In our examples we express the complex 
phase in the CKM matrix in terms of the Higgs fields. 

1.2 Branes at singularities with a broad brush 

Placing D3 branes at singular points in the bulk geometry can reduce the amount of 
supersymmetry and generate chiral matter, providing a rich arena for particle physics 
model building. Models have been constructed on orbifold singularities such as C 3 /Z 3 [7], 
as well as toric and non-toric singularities [25,26j. In this section we provide an introduc- 
tion to recent methods in extracting the gauge theories which arise from branes probing 
toric singularities (much more detailed reviews include [4Tjl4~2"]). 

We focus on conical singularities. This means that we are interested in a Calabi-Yau 
threefold Y whose metric can be locally written in a cone-like form 

ds 2 = dr 2 + r 2 gij dx % daP , (1) 

where < r < 00 and g^ is the metric on a five dimensional Sasaki-Einstein manifold 
The point r = desribes the tip of the cone, where the D3 branes are placed. 
This geometry is non-compact, however one should think of this infinite cone as being 
embedded as a local region within a full Calabi-Yau compactification (for a concrete 
example see [Tll4"8]). 

Generally, an explicit metric for X is not known. However, the gauge theory is 
only determined by the topology of the collapsing hypersurface X, as can be seen in the 

2 As a Calabi-Yau manifold, Y is Kahler and Ricci fiat which implies that X has to be a Sasaki- 
Einstein space [47] . 
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following way: the matter content and superpotential of the D3 branes correspond to the 
matter content and couplings D5 branes encounter when wrapping the two-cycles in the 
non-singular surface X. This is the large volume perspective [4"9lf52| . The gauge theory 
on the D5 branes is determined as in the intersecting brane picture. D5 branes can wrap 
distinct 2-cycles, corresponding to different gauge groups. The intersection of those 2- 
cycles gives rise to chiral matter and triple intersections correspond to superpotential 
terms. 

Focusing on toric Calabi-Yau manifolds allows the powerful techniques of toric geom- 
etry to be used to extract the gauge theory on the branes at the singularity [38,53j (for a 
review of toric geometry see for example [54|). In this context, the six dimensional cone 
can be represented as a T 3 fibration over a convex rational polyhedral cone in R 3 , which 
is parametrised by the normal vectors of each facet. At the facets of the polyhedral 
cone the T 3 fibre degenerates to T 2 , at the edges to S 1 , and at the tip T 3 degenerates 
completely. Imposing the Calabi-Yau condition on Y implies that by an SL(3, Z) trans- 
formation of the torus, all normal vectors can be transformed to Vi = (l,Wi). Being 
rational vectors the basis Wi can then be drawn in the integer lattice of R 2 , forming a 
convex polygon. This is the toric diagram or Newton polytope. 




+ + + + + 



Figure 1: A five- faceted polyhedral cone in K 3 . The normal vectors determine the toric diagram. At 
the facets of the polyhedral cone the T 3 fibre degenerates to T 2 , at the edges to S 1 , and at the tip T 3 
degenerates completely. 

The polygon encodes all the topological information relevant for our purposes. To extract 
particle physics models from a singularity, the analysis always starts with the toric 
diagram of Y. For example, the number of gauge groups is given by twice the area of 
the toric diagram. 

In Section 5 when exploring flavour physics we focus on models which are based on 
toric del Pezzo surfaces. Recall that the del Pezzo surfaces dP n are defined as blow-ups 
at n points (0 < n < 8) of the compact projective space P 2 , the first four of which admit 
a toric description. The associated two dimensional base spaces can be obtained by 
modifying the toric base of P 2 , which is a triangle. For the first toric del Pezzo surfaces 
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the blowing up corresponds to replacing the vertices of the triangle with line segments. 
This is shown for the first del Pezzo surface dP\ in Figure EJ In this case the Calabi-Yau 
3- folds Y are complex cones over the complex two surfaces dP n . 



Figure 2: Left: P2 can be described as a T 2 fibration over a triangle. At the edges of the triangle the 
torus T 2 degenerates to a circle S 1 and at a vertex the fibre degenerates to a point. Right: Blowing 
up a vertex in P2 and replacing it with a line generates the first del Pezzo surface. 

One of the best known examples of a gauge theory from a toric singularity is that 
residing on D3 branes probing the conifold [47J. In this example the five dimensional 
Sasaki-Einstein manifold is T^ l,1 \ and the associated gauge theory is the following J\f = 1 
superconformal quiver gauge theory: 



The ranks of the two gauge groups are equal, the fields A\ 2 transform as (Ni,N 2 ) and 
B l 21 as (N 2 , Ni), where iVj refers to gauge group i. 

From the point of view of model building, a theory with equal ranks is not de- 
sirable. However, by adding fractional branes, the ranks of the gauge groups change. 
In the example of the conifold, adding M fractional branes changes the gauge group 
SU(N) x SU(N) -> SU(N + M) x SU{N). The matter content and superpotential re- 
main unchanged but the theory is no longer conformal [55]. The addition of fractional 
branes allows us to choose the rank of the gauge groups freely and hence to construct 
realistic gauge group content. In addition, the loss of conformality makes realistic model 
building possible. 

Quiver gauge theories 

It can be shown that gauge theories associated to toric singularities are always quiver 
gauge theories |36[l37]. Quiver gauge theories are gauge theories whose matter content 
can be described by directed graphs. Nodes in the graph correspond to gauge groups 
and arrows between two nodes to bi-fundamental matter. Multiple arrows correspond to 






Figure 3: The superpotential and quiver gauge theory of the conifold. 
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multiple copies of that type of matter, i.e. families. The language of quivers is also useful 
in the context of model building. For instance, Higgsing of fields can be implemented 
by collapsing arrows between nodes. The cancellation of non-abelian gauge anomalies 
can be checked by verifying whether the number of incoming arrows and the number of 
outgoing arrows are equal, where each arrow is weighted with the size of the gauge group 
from which it originates or ends. 

The superpotential cannot be read off from the quiver, only all gauge invariant oper- 
ators, which correspond to cycles in the quiver diagram^ Obtaining the superpotential 
leads us to a discussion of dimer techniques. 

Dimer diagrams 

Dimer diagrams encode the complete gauge theory and provide the relation to the toric 
diagram |36tl4"2] . A dimer (also known as a brane tiling) is a graph on a torus T 2 , drawn 
as a parallelogram in M. 2 with opposite edges identified. Each distinct face of the dimer 
corresponds to a gauge group. Edges between faces correspond to bi-fundamental matter 
charged under the associated gauge groups. 

For an edge which divides gauge group faces A and B, the transformation properties 
(i.e. (A,B) or (B,A)) of the bi-fundamental matter are determined by zigzag paths. 
Zigzag paths are directed paths in the dimer which have the property that they turn 
maximally at each node, and they cross once the edges along which they run [391H0] . 
The intersection of two zigzag paths along an edge corresponds to bi-fundamental matter 
whose gauge group transformation properties are determined by the orientation of the 
crossing (cf. Figure H]). The number of common edges between two gauge group faces 
gives the multiplicity of the bi-fundamental matter charged under these two gauge groups. 



Fi gure 4: This crossing of zigzag paths corresponds to bi-fundamental matter transforming as (B,A). 



Zigzag paths form closed paths on the torus. This property is crucial in providing a 
connection with the toric diagram. In particular, the winding numbers of the zigzag 
paths are the inverse slopes of the edges of the toric diagram. 

Vertices in the dimer correspond to terms in the superpotential. All vertices have 
an intrinsic orientation provided by the circulation of zigzag paths around them. Neigh- 
bouring vertices always have opposite orientation, therefore dimers are bipartite graphs. 

3 How to obtain the number of gauge invariant operators is discussed in Appendix iDl 
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The bi-fundamental matter going into a vertex forms a superpotential term with the sign 
of the term corresponding to the sense of orientation, shown in Figure [5] below. 



Figure 5: A part of a dimer showing zigzag paths around two nodes, and the corresponding superpo- 
tential terms. 

With this information in hand, one can obtain both the quiver diagram and superpo- 
tential from the dimer. This is explicitly illustrated for the cases of dPo and dP\ in 
Figure ll.2[ which shows the toric diagram, dual web diagram, dimer, zigzag paths and 
quiver diagram for each of these theories. 




Superpotential for these two vertices: 

W = X13X32X21 — X14X43X32X21 
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However, it should be noted that there can be many different gauge theories which 
correspond to a given singularity, all of which are connected by Seiberg duality [55| I57]. 
The class of theories where the ranks of gauge groups are equal are known as the toric 
phases of the singularity. There can be more than one toric phase associated to any 
given singularity; we will come back to this feature in the context of constraining the 
number of families in Section 3. 

2 Algorithmic view on gauge theories of toric singu- 
larities 

The first algorithmic way to construct the gauge of a toric singularity was presented 
in [53], and the dimer interpretation was introduced in [3B]. Recently, Gulotta [13] 
presented an extremely efficient inverse algorithm to construct the dimer directly which 
allows to construct the gauge theory and provides a visual perspective on many properties 
of the gauge theory. In this section we provide a review of this algorithm, as we shall 
make heavy use of it in the following sections. For proof and further details we refer the 
reader to the original paper [43] . 

The starting point for Gulotta's algorithm is the toric diagram of the singularity of 
interest. When rotated by 90 degrees, the slopes of the edges of the rotated toric diagram 
are equal to the slopes of the dual web diagram. One then embeds the rotated toric 
diagram into the minimal rectangular toric diagram into which it fits. This rectangular 
toric diagram corresponds geometrically to an orbifold of the conifold C. This process is 
illustrated for the case of the third del Pezzo surface <iP 3 in Figure [6j The left hand side 
of this figure shows the toric diagram for dP^, and the right hand side shows the rotated 
toric diagram embedded into a minimal rectangular grid. The philosophy of Gulotta's 
algorithm is to obtain the dimer and gauge theory of interest by partial resolution of this 
orbifold of the conifold. 




Figure 6: LHS: The toric diagram for the third del Pezzo surface dP 3 . RHS: The toric diagram 
rotated by 90 degrees and embedded into its minimal rectangular diagram. 



Since the rectangular toric diagram plays a central role in the algorithm, we would 
like to describe the associated dimer and gauge theory in some detail. An n x m grid 
corresponds to the Z n x Z m orbifold of the conifold C/(Z„ x Z m ). The dimer associated 
to the n x m rectangular diagram consists of n (1,0) and (—1,0) horizontal paths, and 
m (0, 1) and (0, —1) vertical zigzag paths, where the zigzag path (a, b) wraps the torus a 
times horizontally (left to right) and b times vertically (bottom to top). The zigzag paths 
are drawn in alternating fashion such that no paths with the same winding numbers are 
next to each other, e.g. alternating between (1,0) and (—1,0) for the horizontal paths. 
In this dimer there are two types of faces: one corresponding to gauge groups and the 
other to superpotential terms. They are distinguished by the fact that the zigzag paths 
surround superpotential faces in a clockwise or anti-clockwise fashion, whereas for the 
gauge group faces the zigzag paths clash. The intersection of two zigzag paths corre- 
sponds to bi-fundamental matter charged under the two gauge groups at the intersection, 
as in Figure HJ Finally, the different orientations in the superpotential faces correspond 
to different signs in the superpotential terms. Figure [7] shows the toric diagram, dimer 
and quiver associated with a 2 x 3 rectangular grid to illustrate these points. 




Figure 7: Left: Rectangular toric diagram. Middle: The associated chess-board dimer, gauge groups 
are labelled by numbers, superpotential terms correspond to faces that can be surrounded by cycles (the 
orientation determines the sign of the superpotential term). Right: The associated quiver. 

Starting from this rectangular grid, Gulotta's algorithm prescribes how to transform the 
dimer in a way which corresponds to collapsing cycles in the toric singularity, and cutting 
the toric diagram to the desired shape. In this framework, collapsing cycles corresponds 
to merging zigzag paths in an appropriate fashion. Cutting the toric diagram with a 
line of a given slope corresponds to making a zigzag path of precisely this slope. For 
example if we cut off the top left corner of the 2x3 rectangular toric diagram above 
with a line of slope 1, this corresponds to merging a (0, 1) and (1, 0) path in the dimer 
to give a (1, 1) zigzag path. We give precise details on how to merge zigzag paths in 
the next section, and now describe the algorithm to obtain the desired toric diagram by 
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repeated cuttings. 

1. We start with making all paths of slope +1. If the edge removed from the toric 
diagram is an (±n, ±n) vector on the lattice, then n (1, 1) or (—1, —1) zigzag paths 
should be created. 

2. We then cut the toric diagram further with lines of lower and higher slopes, e.g. 
1/2 or 2. This corresponds to creating zigzag paths of higher winding numbers 
such as (1,2). Paths of higher winding number are created by combining already 
existing paths. The order in which these paths are made is given by the Farey tree, 
which is shown in Figure [HJ 

(4,3) 

oaf 

/ 15,3) 

(1,1>H».0»<2,1) 



Figure 8: The Farey tree tells the order in which zigzag paths are made. The part of the tree displayed 
shows the paths with slope < 1, reversing the winding numbers gives the paths with slope > 1. 

3. After cutting the toric diagram with all positive slopes, we cut the edges with 
negative slopes in the same fashion. 

This procedure is continued until the desired toric diagram and dimer are obtained. We 
now return to the issue of merging zigzag paths. 



2.1 Merging zigzag paths 

Let us start with the example of creating a (1, 1) path from the 2x2 grid. Since merging 
paths is a local operation, one can think of this grid as being embedded into any larger 
grid without any loss of generality. 
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Figure 9: Merging a (1,0) and a (0, 1) zigzag path creates a (1, 1) path. 

Combining the (1, 0) and (0, 1) paths gives the desired (1, 1) path, as shown in Figure |HJ 
In a similar manner one can create paths with winding number (1,-1), (—1,1) and 
(—1,-1). Paths with higher winding numbers can be obtained by combining previously 
created paths with further paths of winding number (0, 1) or (1, 0). 

Given two zigzag paths a fundamental topological quantity is their oriented in- 
tersection number. For two paths of winding numbers (m a ,n a ) and (m&, n&) this is 
I a b = iTiaiT-b — n a rrib. On the other hand the total number of intersections is not topo- 
logical, and can change by local deformations of the paths. Note that the number of 
intersections has to be at least as large as the number of oriented intersections. If the 
number of unsigned intersections is greater than the oriented intersection number, we 
say that the two paths have additional crossings. Gulotta's algorithm requires that while 
merging paths we do not generate any additional crossings. Thus while merging paths 
one always performs local operations so as to get rid of any additional crossings. 

As an example, we show in Figure [TU1 additional crossings that can arise while merging 
(1, 0), (0, 1), (-1, 0) and (0, -1) paths to form a (1, 1) and a (-1, -1) path. 
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Figure 10: The additional crossing corresponding to a mass term when a (1,1) and a (—1,-1) path 
are created after each other. 

In this case, to avoid additional crossings one constructs the (1, 1) and (—1, —1) paths 
simultaneously as shown in Figure [TTJ 
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Figure 11: Creating a pair of (1, 1) and (—1, —1) paths simultaneously to avoid additional crossings. 

This crossing has the following interpretation. Recall that the dimer represents the 
massless spectrum of the theory. The circled region in Figure [10] gives a mass term in 
the superpotential. Integrating out the massive fields gives the dimer shown on the right 
hand side of Figure [TT1 

However, not all additional crossings have the interpretation as mass terms in the 
superpotential. The presence of such additional crossings can render the dimer inconsis- 
tent. For example, in the case of multiple paths of one type where I a f, = 0, it is essential 
to avoid additional crossings since otherwise the number of gauge groups does not cor- 
respond to twice the area of the toric diagram. Whether or not a dimer with additional 
crossings is consistent needs to be decided on a case by case For these reasons, 

Gulotta's algorithm requires the removal of all additional crossings, and provides three 
basic operations for merging paths which guarantee the absence of such crossings. 

• Operation I involves making n (±1,±1) paths. One first creates the paths and 
then moves them apart from each other to avoid additional crossings. Some fur- 
ther discussion on removing additional crossings in Operation I can be found in 
Appendix [A] 

• Operation II creates n pairs of (1, 1) and (—1, —1) paths. We saw in Figure [Til 
how to make one pair of such paths without additional crossings. Multiple pairs 
can be achieved by creating one pair and then placing the other pairs parallel to 
the one of the already created paths. An example of this is shown in Figure [T21 
where the dimer of the toric diagram is obtained by creating two pairs of (1,1) and 
(—1, —1) paths. 

4 Examples of consistent dimers with additional crossings are the second toric phase of the zeroth 
Hirzebruch surface Fo and phases 1 and 3 of the third del Pezzo surface dP^,, which are discussed in 
Section 3.3. 
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Figure 12: Example of creating two pairs of (1, 1) and (— 1, —1) paths via Operation 2. The red pair 
is drawn first. Then we draw the second pair (shown in blue) parallel to the (—1,-1) path of the first 
pair. 



Operation III involves creating different numbers of (1,1) and (—1,-1) paths, 
n and m, say. Suppose n — m > 0, then n — m (1,1) paths are created with 
Operation I and the remaining pairs of (1, 1) and (—1, —1) paths are then drawn 
as in Operation II parallel to one of the (1, 1) paths. If n — m < 0, one creates 
m — n (—1, —1) paths and then proceeds as before. An example of Operation III 
is shown in Figure | 




Figure 13: Creating two (1, 1) paths and one (— 1, —1) path by Operation III. We firstly create one 
(1, 1) path (blue) and then create a pair of (1, 1) and (— 1, —1) paths parallel to the blue (1, 1) path. 



The one-to-one relationship between winding numbers of zigzag paths and slopes in 
the web-diagram is essential |39| . Violating this relationship will lead to a dimer that has 
no interpretation in terms of the geometry. This has an important consequence if one is 
trying to construct dimers by merging zigzag paths in a fashion which does not adhere to 
the prescriptions of Gulotta's algorithm. While merging zigzag paths one must not create 
paths which correspond to moving backwards in or out of the Farey tree. One implication 
of this is that zigzag paths are never self-intersecting. For example, combining a (2, 1) 
and a (—1, 0) path corresponds to moving backwards in the Farey tree, while combining 
a (2, 1) and (0, 1) path moves out of the Farey tree. These features have been discussed 
in the context of F in [39J . 
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2.2 Application to toric del Pezzo surfaces 



To illustrate the operations of Gulotta's algorithm we discuss the familiar examples of 
the first three del Pezzo surfaces. Here we are not strict about the order in which the 
paths are created, however it is easily checked that the results do not change, when 
following the algorithm strictly. 

We can embed the toric diagrams of the first four del Pezzo surfaces into the 2x2 
square, as shown in Figure [6] for dP 3 . The dimer associated to this square is the chess- 
board with two paths of each type: (0, 1), (0, —1), (1,0), and (—1, 0). The dimer for dP 3 
is obtained by creating a pair of (1, —1) and (—1, 1) paths, where we use Operation II. 
This process is shown in Figure [T4] Removing additional crossings, as done in Operation 
II, decreases the number of terms in the superpotential. 
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Figure 14: Left: The conifold dimer with the paths to be merged shown in blue and red. Middle: The 
toric diagram for dP%. Right: The dP 3 dimer with the (1, —1) path in blue and the (—1, 1) path in red. 



The superpotential can be read off to be: 
W dP3 ■ 
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Note that every field appears in exactly two terms in the superpotential, once with a 
positive sign and once with a negative sign. In fact, this is a general feature of toric 
singularities, and follows from the bipartiteness of the dimer [37]. We also write the 
superpotential in a matrix form, which we will use in Sections 4 and 5. The matrix will 
play the role of a Yukawa matrix, and the fields in the column vectors will be quarks. 

In order to obtain dP 2 we join the remaining (1, 0) and (0, 1) path from dP 3 , as shown 
in Figure [15j 



17 




Figure 15: Left: The dPz dimer with the paths we will merge highlighted in green. Middle: The toric 
diagram for dP2. Right: The dP2 dimer with the newly created (1,1) path in green. 

Note that the dPi superpotential can be obtained from that of dP% by vevving and 
integrating out the field P53 which was a bi-fundamental under gauge groups 5 and 3. 
These groups were combined in the merging of the zigzag paths. This is a general feature 
when cycles are collapsed. 
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dPi is obtained by collapsing another cycle, this involves generating a (—2, 1) path. We 
note that even when paths of higher winding number are created, the structure of the 
dimer only changes locally. 
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Figure 16: Left: The dP 2 dimer with a (—1,0) path highlighted in red. Middle: The toric diagram 
for dP\. Right: The dP\ dimer with a (—2, 1) zigzag path (red). 
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Generating a (1, —2) path out of the existing (1, —1) and (0, —1) path leads us to the 
dimer of dP . 



1 




.4- 


1 


+ 


3 , 


+ 


2 1 


+ 


3 , 


y " 


1 


's « 


- 3 


+ 


t 2 




1 1 


+ 


1 








1 




1 




• -. 


1 


+ 


3 1 


+ 


2 1 


+ 


3 


y ' 


1 


's « 


3 


+ 


2 




1 


. + 


1 


-VI 




1 



Fi gure 17: Left: The dP\ dimer. Middle: The toric diagram for dP$. Right: The dPg dimer with the 
(1, —2) path shown in blue. 

The superpotential for dP exhibits an SU(3) flavour symmetry under the rotation of 
families. 
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Connecting Gulotta's dimer with the traditional dimer 

While the zigzag path diagram introduced by Gulotta contains all the information en- 
coded in dimers, their visual appearence is quite different. Here we discuss how one 
can reassemble dimers from dimers constructed with Gulotta's algorithm. This process 
is shown in Figure [THJ Beginning with the zigzag path diagram corresponding to the 
conifold C, we label the gauge groups and superpotential terms. We then draw lines 
between the centres of the superpotential faces. The zigzag paths can then be bent, to 
make the structure resemble the dimer diagrams of |37j . 
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Figure 18: The connection between zigzag path diagrams and dimers. 

2.3 D7 branes 

In local models D7 branes enter by wrapping non-compact holomorphic four cycles which 
pass through the singular locus. In a given singularity the number of different four-cycles 
is given by the number of external points in the toric diagram. In a compact model the 
D7 branes would wrap cycles which would extend into the bulk region of the Calabi-Yau 
and have non-vanishing gauge coupling. The presence of the D7 branes leads to three 
new kinds of string states, 37, 73 and 77 strings which we discuss in turn. 

For every bi-fundamental 33 field there exists a pair of 37 and 73 states which interact 
with the 33 state via a 33-37-73 coupling. Thus, when families of 33 states exist, there 
can be multiple 37 and 73 states with which they interact. Dimer techniques for the 
computation of the 37 and 73 spectrum and superpotential have been developed in [58J. 

If there are multiple D7 branes there will be D7-D7' strings leading to 77' matter 
fields which interact with 37 fields via a 37-77'-7'3 coupling in the superpotential. Un- 
fortunately, since the 77' fields are supported on noncompact cycles dimer techniques 
are not sensitive to the full details of their spectrum and interactions. 

For general geometries, the 33-37-73 couplings are diagonal in flavour space. If the 
singularity possesses isometries there is no longer a unique correspondence between one 
pair of 37 and 73 states and 33 states. For instance, in the case of the C 3 /Z 3 orbifold [7] 
the orbifold action is Zj — > Zje^f . This implies that any equation ^ . (3jZj = is invariant 
under the orbifold action, leading to a freedom in the choice of holomorphic four cycles 
that D7 branes wrap. This leads to a more general form of the superpotential, in which 
a given 33 state might couple to an arbitrary 37 state with the strength of the coupling 
being determined by /3j. This freedom arises from the SU(3) global symmetry associated 
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with being able to perform SU(3) transformations on the coordinates Z{. Moving away 
from orbifolds, these arguments can be generalised [58J to show that the presence of a 
continuous global symmetry implies off-diagonal interactions between 37 and 33 statesj§ 
In dPo the singularity has an 577(3) isometry, in dP\ there is an SU(2) isometry and for 
the higher del Pezzo surfaces there are no non-abelian isometries. 

In the context of models which involve 37 states [7J|S] this implies that at the dP® 
singularity the Yukawa matrix for the down generation is a generic 3x3 complex matrix, 
while for dP\ one has mixing between two of the three generations. For the higher del 
Pezzo surfaces the Yukawa matrix is diagonal. These facts shall play an important role 
in our discussion of flavour physics in Section 

3 Families from toric singularities 

Families are matter in gauge theories with the same quantum numbers, which in quiver 
gauge theories corresponds to the number of arrows between two gauge group nodes. In 
this section we explore the bounds on the number of families which exist in gauge theories 
probing toric singularities. We first discuss families in the context of Gulotta's algorithm 
and find that there is an upper bound of three families. After that, we include the 
possibility of addtional crossings and see that there is a unique exception to this bound, 
namely the zeroth Hirzebruch surface, one phase of which has four families. Finally we 
discuss the effects of Seiberg duality and how the bound may be circumvented. 

3.1 Families a la Gulotta 

In this section we consider the effects which the operations of Gulotta's algorithm have 
on the number of matter fields and their couplings. We will show that no more than six 
matter fields exist charged under any given gauge group. Since the toricity condition 
implies that three of these are fundamental and three anti-fundamental, there can be no 
more than three families. 

The starting point for obtaining the gauge theory associated with a given toric singu- 
larity is the smallest rectangular diagram in which the toric diagram can be embedded, 
corresponding geometrically to an orbifold of the conifold. In such a theory there is at 
most one arrow connecting any two gauge groups; every second square in the zigzag 
path diagram corresponds to a different gauge group, and there is not more than one 
field between any two gauge groups. An example is shown in Figure [7] which shows a 
2x3 rectangular toric diagram, its zigzag diagram and its quiver. Our starting point 

5 We thank A. Uranga for pointing this out. 
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will therefore always be a quiver with only one family§| 

While the effects of Gulotta's algorithm are local in the dimer for the gauge group 
faces and couplings, they are global from the point of view of zigzag paths. However, 
the changes in the global nature of the zigzag paths only affect the gauge groups and 
couplings in the region where zigzag paths are being merged. Thus, while we present our 
results based on a small part of a larger dimer, the embedding this into a larger dimer 
will have no effects on our results. We now proceed systematically through the effects of 
the operations of Gulotta's algorithm. 

• Operation I creates n (1,1) paths. Consider the creation of a single (1, 1) path, 
which corresponds to joining two gauge groups and changing two quartic couplings 
to cubic couplings, as shown in Figure [THJ 
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Figure 19: The effect of joining a (1,0) and (0,1) path (Operation I) on a part of a larger dimer. Gauge 
groups 4 and 5 are joined, and two quartic couplings become two cubic couplings. 

When creating multiple (1, 1) paths additional crossings are inevitable, and have 
to be removed. While this leads to global changes in the structure of the zigzag 
diagram, locally the situation always resembles the figure above since removing 
the additional crossings creates a situation locally the same as in the above figure. 
This is illustrated in Figure [32] in Appendix [A] From Figure [19] it is clear that 
there are 6 matter fields charged under gauge group 4. The toric condition implies 
that three of these transform in the fundamental representation of gauge group 4, 
and the other three as anti-fundamentals. There can therefore be a maximum of 
three families arising from this operation. 

• Operation II creates n pairs of oppositely winding paths. As an example in 
Figure [20] we present the case of a pair of (1,-1) and (—1,1) paths. By this 



3 This does not apply to the conifold itself, which has two families. 
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procedure two pairs of gauge groups are joined, 6 with 12 and 4 with 11, and five 
quartic couplings become two cubic and one sextic coupling. 




Figure 20: The effect of Operation II, in this case creating a (—1, 1) and a (1, —1) path. Two pairs of 
gauge groups are joined, and five quartic couplings become two cubic and one sextic coupling. 

By the same logic as above, not more than two families can be created here. If 
we wish to create n pairs of oppositely oriented paths we may do so according to 
the prescription in Section [2j From Figure [12] in that section it is clear that, in 
creating further zigzag paths the number of fields charged under any given gauge 
group does not change. Since every gauge group has four matter fields charged 
under it, there will still be no more than two families. 

• Operation III creates different numbers of (1, 1) and (—1, —1) paths, m and n, 
say, where m > n. The first step is to create m — n (1, 1) paths. This is just 
Operation I, and creates no more than 3 families. One then removes n pairs of 
(1,0) and (—1,0) paths, and n pairs of (0,1) and (0,-1) paths. This does not 
change the number of families. Finally, we create n pairs of adjacent (1, 1) and 
(—1,-1) paths. To avoid additional crossings, these should be created parallel to 
one of the already existing (1, 1) paths from the first step. By creating the paths in 
this manner, no extra families are created, and so for Operation III the maximum 
number of families obtainable is three. 

Gulotta's algorithm involves repeated applications of these operations. For example 
zigzag paths with higher slopes are generated by applying multiple Operations I, II or 
respectively III. What are the effects of this? The only non-trivial effect appears when 
such additional operations affect a local region that has already been modified, as this 
is the only way to modify one gauge group that has been modified already. Typically, 
joining gauge groups together generates additional crossings which after removal do not 
lead to extra families. 
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For example, consider the zigzag diagram on the left of Figure |5T] which has a single 
(1, 1) path, created by Operation I. If we tried to add a (—1, —1) path to this diagram to 
form a longer gauge group "tube" using the same operation, we would obtain additional 
crossings in our diagram, as in the right hand side of Figure [2TJ Gulotta's algorithm 
prescribes the removal of these additional crossings in such a way that no further families 
are introduced. 
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Figure 21: Left: A zigzag diagram with a single (1, 1) path and two highlighted (—1,0) and (0, — 1) 
paths which will be merged. Right: The same diagram after the paths have been joined. There is one 
(1, 1) and one (— 1, —1) path, but with additional crossings. 

The only case where a non-trivial modification of the gauge group faces takes place 
without introducing additional crossings is by applying Operation II or III after Opera- 
tion I. An example of this is in the transition from dP 2 to dP\ as shown Figure [T6j Doing 
this generates at most 3 families. 



3.2 Additional crossings 

While the dimers obtained with Gulotta's algorithm do not include additional crossings, 
their presence does not necessarily render the dimer inconsistent. In this section we 
therefore consider the possibility that the paths in the zigzag diagram have additional 
crossings. We will show that the maximum number of families is four and that there is 
only a single toric gauge theory which has four families: the zeroth Hirzebruch surface 
F . 

Once we allow the presence of additional crossings, there is no reason for arbitrarily 
long gauge group tubes not to exist, such as shown in the right hand side of Figure I2T1 
However, we may only create straight gauge group tubes, not ones that 'bend'. Bending 
gauge group tubes lead to inconsistencies in the dimer, as discussed in [39] and mentioned 
in Section 2.1. An example of such a path is a (1, 1) path made from the combination of a 
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(2, 1) and a (—1, 0) path. The presence of such paths violates the correspondence between 
zigzag paths and slopes in the web diagram, and the associated dimer is inconsistent. 

Since the gauge group tube cannot have addtional branches it must have only one 
beginning and one end. There are then two possibilities: that the gauge group tubes in 
which we are interested are orthogonal to each other or they are parallel. The first case is 
heuristically shown on the left of Figure [22] where the orange and purple tubes represent 
two different gauge groups. The arrows show the matter fields which are bi-fundamentals 
of the orange and purple gauge groups. It is clear that a maximum of four families is 
possible in this case. While the diagram shows a square chess-board, it is possible to 
generalise this to a rectangular one. 

The right hand side of Figure [22] shows the case of parallel gauge group tubes. In 
this case there are superpotential mass terms for the bi-fundamental matter between the 
orange and purple gauge groups which must be integrated out. After this there will be 
at most one family of chiral matter, this occurs if the common length of the gauge group 
tubes is odd. 




Figure 22: Left: A zigzag diagram demonstrating the presence of four families with orthogonal gauge 
group tubes. The purple and orange tubes are different gauge groups, with arrows highlighting the 
bi-fundamental matter between purple and orange gauge groups. The green arrows show a zigzag path 
with an inconsistent winding that would be present in this dimer. Right: A zigzag diagram with parallel 
purple and orange gauge group tubes. The matter fields are non-chiral. 

The schematic four family case from the left hand side of Figure [22] is not easy to realize 
since creating orthogonal tubes leads to inconsistent paths, as mentioned in the previous 
section. The green path in the above figure has winding (0, 0), which is inconsistent with 
the correspondence to the web diagram. The only situation that does not contain such 
inconsistent paths is the most simple one, which is the zeroth Hirzebruch surface Fo, as 
illustrated in Figure [23] 
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Figure 23: From Left to Right: Toric diagram of the zeroth Hirzebruch surface. The dimer obtained 
with Gulotta's algorithm, containing two families. Four family Fo dimer, obtained by allowing for 
additional crossings. A schematic zigzag diagram highlighting how the four families arise. 

We stress again that the 4-family scenario is only possible by allowing for additional 
crossings in the zigzag diagram. Those additional crossings can be avoided while con- 
structing the dimer according to Gulotta's algorithm. Since we have shown above that 
Gulotta's algorithm leads to no more than three families, this implies that any singu- 
larity which has a toric phase with four families must also have one with less than four. 
This leads us to some comments on Seiberg duality. 



3.3 Seiberg duality 

There are various gauge theory descriptions for branes probing a toric singularity, all of 
which are connected by Seiberg duality [56"l l5? l l59 ] . Gulotta's algorithm provides us with 
one gauge theory probing a given singularity in a toric phase. There are two ways of 
obtaining different toric phases: 

1. Start with a gauge theory which corresponds to a toric phase of the singularity. 
Then Seiberg dualise one of the gauge groups such that the gauge theory remains 
in a toric phase. 

2. Allow for consistent additional crossings in the construction of the gauge theory, 
as was done previously for Fo. 

We wish to consider Seiberg dualising a gauge group S of a toric gauge theory subject 
to the requirement that the dual theory also be in a toric phase. Say that the rank of 
the gauge groups is N c and that it has Nf flavours. The dualised gauge group has rank 
N' c = Nf — N c after the duality. Since toricity requires N' c = N c = N, the number of 
flavours must be 2A" and hence there are two fields transforming in the fundamental and 
two in the antifundamental representations of S. 

In the quiver gauge theory, dualising one particular node replaces all outgoing fields 
Qsa, which transform as (Ng, N a ), and incoming fields Qbs, which transform as (Nf,, N$), 
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by their dual fields q a g (N a , Ng) and q$b (Ng, iVj,). This implies the arrows corresponding 
to those fields in the quiver are reversed. In the dual theory one also adds all possible 
mesons made out of bound states of the original fields, My = QisQsj- This can easily be 
visualised in the quiver, as shown in the left hand side of Figure [2H The superpotential 
is modified to 



where all indices are suppressed. 

In the dimer we can form a picture of Seiberg duality following [57J . The right hand 
side of Figure [241 shows a local region of some larger dimer with a gauge group face S with 
2N flavours and the effects of dualisation. Dualising changes the white nodes around 
gauge group S to black and black to white, which reverses the transformation properties 
of the bi-fundamental matter. The presence of Seiberg dual mesons is encoded in extra 
edges between gauge groups 1 to 4 and gauge group S which are shown in red. 



Figure 24: Seiberg Duality on Quiver. Left: A local part of a larger quiver. Right: The quiver 
after dualising node S. The fields are reversed (blue) and mesons are added (red). Seiberg Duality 
on Dimer. Left: A local part of a larger dimer. Right: The same face after Seiberg dualisation. The 
white and black nodes are reversed and there are new mesons shown in red. Black and white nodes 
correspond to opposite signs in the superpotential. 

We can now write down the effect of Seiberg duality on zigzag diagrams. The left 
hand side of Figure [25] shows a local section of a zigzag diagram which corresponds to 
the local section of the dimer shown in Figure EH To dualise this zigzag diagram we 
must introduce Seiberg dual mesons. In order to reverse the (Ng, A?j) fields we also need 
to change the orientation of the zigzag paths surrounding face S. We can do this by 
twisting the zigzag paths. The net effect is shown in the right hand figure of Figure [25] 
and is local in the zigzag diagram. 



W duiiy (qq, M, $) = W OTig (M, $) + Mqq , 



(6) 
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Figure 25: Le/t: A local part of a zigzag diagram which will be Seiberg dualised. Right: The effect of 
dualising the left hand zigzag diagram. The orientation of the paths around face S changes and twists 
in the zigzag paths are introduced. 

Another way of obtaining toric phases of the gauge theory is by allowing for additional 
crossings of zigzag paths. Specifically, these are additional crossings that are not of mass 
type and not between paths of same winding number. The canonical example of this is 
the creation of a pair of paths with opposite winding; unlike in Operation II they are not 
created at the same time but via Operation I performed twice. The crucial point is that 
the two paths are created locally distinctly such that they do not generate a mass term 
unlike in Figure [TUl As an example, in Figure [261 we show phase 3 of dP 3 [57| . which is 
Seiberg dual to the dP^ gauge theory obtained in Section 12.21 from Gulotta's algorithm. 
Phase 3 is obtained if one creates the (—1, 1) and (1, —1) paths of dP% after each other 
and in a locally distinct fashion. 

We should note that using this method only a limited number of toric phases can 
be found, whereas all toric phases may be obtained by using the twisting procedure 
described earlier in this section. 
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Figure 26: Left: Dimer for phase 3 of (IP3, obtained by allowing for additional crossings. Right: Dimcr 
of phase 4 of dP-$, obtained by dualising gauge group "2." 

What are the implications of this discussion for the numbers of families? The toric 
condition limits us to dualising gauge groups with 2N C flavours. This means that there 
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must be two fields transforming in the fundamental and two in the anti-fundamental 
of the gauge group which will be dualised. Therefore, there can be a maximum of four 
Seiberg dual mesons and a maximum of four families from Seiberg dualisation is possible. 
Given the arguments we have presented this section and the previous one, Fo is the only 
example where Seiberg duality leads to four families. 

3.4 Obtaining further families 

There are a number of ways to obtain further families, thereby evading the bounds we 
have derived. Our bound is based on the assumption of toricity. In terms of zigzag paths 
this implies that there are no self-intersecting or inconsistently constructed zigzag paths. 
From a field theory perspective one can Seiberg dualise gauge groups such that one is no 
longer in a toric phase. It is known that the non-toric phases often exhibit more than 
four families [52} 160]. 

There is currently no dimer-based technique available for the investigation of non- 
toric phases. Using zigzag diagrams one can dualise gauge groups by twisting zigzag 
paths. It may be that such twist operations could be used to understand non-toric 
phases, possibly leading to more families. We find the prospect that zigzag diagrams 
could lead to a greater understanding of the non-toric phases to be an exciting one. 

Another way of circumventing the bound is by Higgsing fields in the quiver. With 
Higgsing one can in principle create bendy tubes which could lead to the proliferation of 
families. 

4 Yukawa couplings and masses at toric singularities 

In this section we explore the structure of Yukawa couplings for gauge theories at toric 
singularities. The models that we shall discuss shall have unequal gauge groups at 
the nodes and hence non-conformal. In the absence of supersymmetry breaking they 
would flow to strongly coupled superconformal field theories in the infrared. We note 
that at scales below the SUSY breaking scale the renormalization group flow would be 
significantly modified, as in [8] we would require to arrange for scales such that the gauge 
couplings are weak at the SUSY breaking scale. We show that for models capable of 
accounting for the multiplicity of families, one quark mass is zero at tree level, whereas 
the other two are generically non-zero and can differ hierarchically. 
Let us start with the quark sector of the MSSM superpotential, 

w d + viQWL (7) 

where y u and y<j are 3x3 Yukawa coupling matrices. Although the MSSM has only two 



29 



Higgs fields it is natural in string phenomenology to have more, as for instance in the 
models of 0125] ■ The values of the Yukawa couplings are then determined by the values 
of the Higgs vevs and those of any moduli which may appear in the superpotential, and 
hence depend on the geometric realisation of the model. Toric singularities are special 
in this regard since there are no complex structure moduli fields which appear in the 
superpotential, and the mass hierarchies and Yukawa couplings depend only on the vevs 
of the extended Higgs sector. The superpotential for an MSSM-like model at such a 
singularity is then best written as 

W = Q l L c^H k u y + Q\c t] (R k d )d\ (8) 

where c^- and are functions of the extended Higgs sector. 

As an explicit example, consider the Standard-Model dPo model of [7], shown in 
Figure E3 




Figure 27: The quiver for the Standard Model realised at a dPo singularity, taken from [7J. 
The up-quark Yukawa couplings in this model are determined by the dPo superpotential: 
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where the matrix of Higgs fields is the Yukawa matrix. It is this matrix which determines 
the flavour structure and masses. The squares of quark masses associated with the 
Yukawa matrix are given by the eigenvalues of 



M = Y.YK 



(10) 



For dPo the mass eigenvalues are which appears once and |ATi2| 2 + \ Y\2\ 2 + I-Z12I 2 , which 
appears twice. Hence the dPo singularity lacks sufficient structure to provide realistic 
Yukawa couplings. As a resolution, this model was embedded into dP\ [8J and the 
eigenvalues were found to be 
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We find that for the higher del Pezzo surfaces the mass eigenvalues have the same 
structure (0, m, M) where m 7^ M. The zero eigenvalue can be explained as follows. 
Firstly, the toric condition implies that every field appears twice in the superpotential, 
once with a positive sign and once with a negative sign. The determinant of the Yukawa 
matrix is the difference of two terms which are equal, each term being the product of all 
the Higgs fields. Since the matrix M is the product of Y and its Hermitian conjugate, 
the vanishing determinant of Y guarantees the existence of a zero eigenvalue of M. 

In the following section we show that the vanishing determinant is a generic feature 
of toric singularities, subject to physically motivated constraints regarding the choice 
of fields labelled as quarks. We comment further on the possible mass hierarchies in 
Section [OJ 

4.1 Mass matrix for the rectangular grid 

In order to discuss universal features of the mass matrix in toric singularities, we start 
with the father theory of all toric singularities, the rectangular grid. The dimer associated 
to the rectangular grid is a chess-board (cf. Figure [7]). Our strategy is to convert any 
statement true for the grid to a corresponding statement in the daughter theories made 
by partial resolution of the singularity. While resolving the singularity we shall allow for 
generic higgsings, even those which take us to non-toric singularities 

For the rectangular grid all interactions are quartic and each field appears exactly 
twice in the superpotential. Before proceeding we therefore must clarify the meaning 
of 'quarks' and 'Yukawa couplings'. A Yukawa coupling in this context is a term in 
the (quartic) superpotential which involves exactly two quarks, one left-handed and one 
right-handed. We require that there are no terms in the superpotential which involve 
more than two quarkgj (for both left and right handed). 

This implies that the set of the left-handed (and respectively right-handed) quarks is 
part of a perfect matching of the dimer. This leads to an upper bound on the number of 
fields we may designate as 'quarks'. In particular, for annxm chess-board the maximum 
number of 'quarks' we may choose is n • m. 

By Higgsing one can break two gauge groups to the diagonal gauge group, leading to 
family replication from quarks which were originally at different places in the quiver and 
transforming under different gauge groups. An example of this is given in Appendix O 
To take this possibility into account, we adopt a broad definition of quark in this section: 
the quarks are just some specific choice of fields in the chess-board, free of any constraints 
of transformation under the same gauge group. They are subject to the constraint above, 

7 A set of fields with this property has been discussed in studies of formal aspects of dimers, and is 
dubbed a perfect matching. 
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namely that not more than two fields in any of the quartic superpotential terms are 
quarks. 

There is one final consistency condition in choosing the quarks. This is that the quarks 
should be chosen as connected or alternating lines, such as in Figure [28j Choosing the 
quarks via disconnected lines, as in Figure [301 makes it impossible to vev the remaining 
fields in such a way that one ends up with a common gauge symmetry under which all 
quarks transform. 

In what follows we deal systematically with the three distinct possibilities that arise 
from the above discussion. The first is that we choose the maximal number of 'quarks', 
n ■ m. The second possibility is that there is less than the maximal number of quarks. 
Finally and for completeness we discuss the seemingly unphysical case where the quarks 
are chosen in a disconnected manner. 



A: Maximal number of quarks 



To illustrate our argument we consider the specific case of the 2x2 chess-board. Our 
result also holds more generally for the n x m chess-board. Figure [28] shows a choice 
of perfect matchings which corresponds to those chosen for the del Pezzo surfaces in 
Section I2T2"} apart from the changes in quarks X 86 , X 8 i, and X 2 $- In the case of the del 
Pezzos we had the corresponding quark Xsg. 
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Figure 28: Here we display the two perfect matchings which correspond to the choice of quarks in the 
Yukawa matrix Eq. 1121 



The Yukawa matrix in this case is 
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The determinant of the Yukawa matrix vanishes and there is one vanishing eigenvalue. 



B: The case of three quarks 



If we do not choose the maximum number of possible 'quarks' then the situation changes 
slightly. We consider the case of three quarks in the 2x2 chess-board. 
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Figure 29: Choosing three quarks in the grid. 
We can write the superpotential as follows 
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The determinant of the Yukawa matrix is non-zero: 

det(Y) = XyjX^X^iXx^X ^2X^1 — X 67 X 72 X 8e ) 



(13) 



(14) 



One can see by vevving fields (X72) = m\ and (Xgg) = (joining gauge groups 2 & 7 
and 6 & 8 as in dP 3 ) the bracket in the determinant simplifies to X^X^Xgi — m\m2X§2- 
One finds that the superpotential has mass terms for the fields X 62 and X 2 q. Integrating 
them out leads to the vanishing of the determinant. 

We observe that although we started with a non-vanishing determinant in the first 
place, vevving and integrating out lowers the rank of the Yukawa matrix, generating one 
zero eigenvalue. Again, our argument should generalise to larger numbers of quarks on 
larger chess-boards. 



33 



C: Disconnected choice of quarks 

The final possibility is to take a disconnected choice of quarks. It has been mentioned at 
the beginning of this section that in this case there is no way to Higgs the theory such 
that the quarks share a common gauge symmetry. We speculate it may be possible to 
construct realistic theories based upon orientifolding. 
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Figure 30: A disconnected choice of quarks. 
For the disconnected choice of quarks in Figure [30] the Yukawa matrix is 

W = X14X31X45X53 + X23X36X42X64 — X17X23X31X72 — X14X28X42X81 

+X17X58X75X81 — X53X67X75X81 — X45X58X64X86 + X28X67X72X86 
/ X23 \ / — X 17 X 72 X 4 2 X 64 \ / X31 > 

_ x 53 X 14 X 45 — X 67 X75 x 36 
X58 X17X75 — X45X64 Xsi 

V X28 / \ — X14X42 X67X72 / \ Xs6 / 

In this case the determinant of the Yukawa matrix does not vanish: 

det Y = X^X1 2 X1 5 Xq 4 — X14X17X42X45X64X67X72X75 . (16) 

There are many different ways to embed non-maximal numbers of quarks once the 
requirement of connectedness is dropped. However, given a non-maximal non-connected 
choice of quarks, it is always possible to embed them into a maximal embedding, by 
completing the perfect matching associated with the maximal choice. This makes it 
clear that whatever the non-maximal choice is, the associated Yukawa matrix can be 
found by removing rows and column from the maximal Yukawa matrix. Whether the 
determinant vanishes before Higgsing depends on the choice of quarks, and it is possible 
to find examples with and without zero eigenvalues. 
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4.2 Mass hierarchies 



The hierarchy amongst the quark masses observed in nature is highly puzzling. In 
string phenomenology these masses are tiny compared to the string scale in most models 
and it is a very delicate problem to generate the correct quark masses dynamically. 
Ultimately the size of these masses is related to cycle volumes and the effect of higher 
order corrections needs to be carefully analysed. Thus, zero quark masses for the lightest 
generation seems to be a natural starting point. The zero quark mass which we have 
found to be generic in toric singularities is therefore very attractive. Let us now explore 
whether the other two eigenvalues can be hierarchical, neglecting higher order corrections 
(e.g. from the Kahler potential). 

In the case of one zero eigenvalue the characteristic polynomial of the mass matrix 
simplifies to a quadratic equation. For dP\ this leads to the non-vanishing eigenvalues: 

I $61 " 



Ai 




A 2 



(17) 

For the field theory description to be valid, the vev for $6i has to be (much) smaller 
than the cutoff A. The physically interesting limit of $6i/A <C 1 leads to a hierarchy 
between the two non- vanishing eigenvalues, unless X,Z < F [8]. The expression for the 
mass eigenvalues for the higher toric del Pezzo surfaces are complicated. However, it is 
easy to isolate conditions under which they are hierarchical. In dP 2 , the discriminant of 
the solution of the quadratic equation is given by: 

A 4 |X 12 | 4 (A 2 - |$ 61 | 2 ) 2 + (|F 64 | 2 (A 4 - A 2 |* 42 | 2 ) + |Z 14 | 2 (A 4 - |$ 61 | 2 |* 42 | 2 )) 2 

-2A 2 |X 12 | 2 (A 2 - |$ 61 | 2 ) (|F 64 | 2 (A 4 - A 2 |* 42 | 2 ) + |Z 14 | 2 (-A 4 + |$ 61 | 2 |v& 42 | 2 )) 

> (A 2 |X 12 | 2 (A 2 - |$ 61 | 2 ) - (|F 64 | 2 (A 4 - A 2 |* 42 | 2 ) + |Z 14 | 2 (A 4 - |$6i| 2 |*4 2 | 2 ))) 2 • (" 

In the physical scenario of small vevs for <3> 61 and \I/ 42 , to obtain a hierarchy we also have 
to impose the following condition: 
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This is satisfied if any of the ratios in the denominator are large. Clearly there are more 
directions in the g?P 2 moduli space compared to the dPi moduli space that generate a 
hierarchy. In the case that the discriminant equals its lower bound and satisfies the 
above constraint, the eigenvalues are given by the hierarchichal structure: 



I I-1ZI . |-04| . |-14| , 
!> 2 1 IA„,|2,_. ,„ IVT,.„|2 ,„ IA„ |2|,T,.„|2 , _ ,„ ) 




35 



Similarly, the hierarchy in <iP 3 is at least 

IV- 12 i |V 1 2 i JP53P17 12 
_ . |A 12 | +1^41 +^p-|Z 14 |, 

1,2 ' '*6l| 2 |y 12 , l*42f|V 12 | Igglflggflrf lp |2 1 ' 

-p— |Al2| H JT~ I ^ 64 1 H X s l Zl4 l ' 

Whether this hierarchy generalises to other singularities depends on the particular 
Higgs structure and has to be examined on a case by case basis. However, the eigenvalues 
will generically be different and hence such limits ought to exist. 



5 The CKM matrix 



In this section we explore the flavour structure of gauge theories at toric singularities, 
with the goal of obtaining models which have a realistic CKM matrix. We begin with 
a short review of flavour mixing in the Standard Model, primarily to set conventions. 
When written in terms of quark fields which are gauge eigenstates, the flavour structure 
is encoded in two 3x3 complex matrices, the Yukawa matrices Y u and Yd of the up 
and down generations. The change of variables to mass eigenstates involves rotations by 
independent unitary transformations for the left and right handed quarks {ul —> V u Ul, 
ur — > V u Ur, d,L — > Vddh and df> — > Vddn). These transformations diagonalise the Yukawa 
matrices (Y u — > V^Y U V U = D u , Yd — > VjYdVd = Dd) and induce flavour changing 
processes via the weak current. The coupling between various generations is given by 
the CKM matrix 

VcKM = K f ^. (22) 
A useful method to compute the CKM matrix is to consider the Hermitian matrices 
Y-vXu an d Yd¥<i i ■ Vu and Vd then correspond to the unitary matrices associated with the 
similarity transformations which diagonalise Y U Y^ and Y^yJ. Experimentally the CKM 
matrix is EH 
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To leading order in each entry of the matrix this can be parametrised hierarchically as 
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0.2. Another useful parametrisation is by three angles and a complex phase, 
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using the convention that q = cos#j and Sj = sin^j. For now we shall focus only on 
obtaining the appropriate hierarchies in the CKM matrix, and take its entries to be real 
which reduces the parametrisation (1231) to the Euler angle decomposition of orthogonal 
matrices. Then, each of the matrices in Eq. [25] induces a mixing between two generations. 

As discussed in Section HI models arising from branes at singularities typically involve 
multiple Higgs fields. We find that the CKM matrix is very closely related to the ratios 
of the vevs of these fields. A study of the complete Higgs potential is beyond the scope 
of this paper, it being a combination of various effects such as D-terms, supersymmetry 
breaking and radiative corrections. For the present analysis we take a more phenomeno- 
logical approach, allowing arbitrary vevs of the Higgs fielda§ and attempting to identify 
the regions in field space which give a realistic CKM matrix. Following the philosophy of 
the bottom-up approach, our analysis can be thought of as trying to answer the question: 
is there enough structure in the singularities to allow the construction of models with 
realistic flavour physics? In our analysis we distinguish between two broad classes: 

1. Models with up and down quarks coming from different sources, i.e. D3-D3 and 
D3-D7 states, and 



2. Models with both up and down quarks coming from D3-D3 states. 



5.1 Up and down quarks from different sources 

We begin with the discussion of the CKM matrix associated to the SM-like model (its 
dP realisation is shown in Figure |27|) . The dPo realisation has the unsatisfying feature 
that the mass of the two heaviest quarks in the up generation are equal. In addition, 
as discussed in Section 2, the SU(3) flavour symmetry of the dPo geometry allows for 
a generic Yukawa matrix in the down generation. This implies the from of Vd is un- 
restricted. As a result the form of the CKM matrix is arbitrary and the appropriate 
hierarchies can be obtained by tuning the vevs of the Higgs fields. Next, let us discuss 
the more realistic dP\ model. 

8 We take care that the vevs of the Higgs fields are always below the cutoff scale as is required for 
the validity of the effective field theory. 
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The dPi model 



The superpotential for the up quarks in the dP\ model was given in equation HI The 
associated unitary matrix for the up quarks V u is 
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X 12 \ 2 + 



\Z 12 \ 2 . One of the features of the matrix is that one of the entries vanishes. We shall 
find that such zero textures are quite common. The basic structure of V u is similar 
to that of dPo, and V u for dPo can be obtained from Eq. [26] by setting $gi = A, but 
the form of Vd is lot more restricted than in the dP model. As discussed in Section 2 
the SU(2) flavour symmetry allows for mixing between two generations of quarks and 
the corresponding matrix Vd induces mixing between only two of the three generations. 
Now consider D7 brane wrappings which correspond to mixing between the medium and 
heavy generations. Explicitly, we take the down quark Yukawa matrix to be 



Y, 



( m d 




The corresponding matrix Vd is 













m s ci 


+ m h s\ 




- m s )s 2 c 2 


(m b - 


- m s )s 2 c 2 


m b c 2 + m s s\ 




( 1 ° 


^ 




v d = 


c 2 








V o s 2 


C 2 J 





\ 



(27) 



which is of the form R 23 . 

Now, let us examine whether it is possible to obtain a realistic structure of the CKM 
matrix. For this purpose we consider V u in the limit j^, J^x 12 ^ ^' * n ^o 
leading order 
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Note that the Hermitian conjugate of this matrix has precisely the same form as the 
product two of the three rotation matrices in 
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with Si « — V A1 1 2 and s 3 rj -f 12 -. 

Thus the CKM matrix Vqkm = Vu^d has the right structure to describe mixing 
between all the generations. To obtain the correct hierarchical flavour mixing, we require 

AY12 Z12 3 2 /oi\ 

~ e , s 2 ~ e . (31) 



X 12 $61 -^12 

Note that in this limit the mass eigenvalues are hierarchical. Therefore dPi is attractive 
from the perspective of both mass hierarchies and flavour mixing. 



Higher dPs 

One can also consider extending the model to the higher del Pezzo singularities. The 
absence of isometries in the geometry implies that at tree-level the mass matrix for the 
down quarks is diagonal, and so the matrix Vj is equal to the identity matrix, implying 
Vckm = Ki- The V u matrices are complicated and we do not present their explicit 
form. One feature is that the texture zero present in dPo and dP\ is absent. It can 
also be shown that there are relationships relating the second-third generation mixing to 
the first-second and first-third which force the second-third generation mixing to be of 
smaller magnitude than the first-third § Thus it is not possible to obtain a fully realistic 
structure at tree level when the model is extended to the higher del Pezzo surfaces. 

However, we note that there is a consistent truncation where the only no n- vanishing 
mixings induced by V u are those between the first and second generations. This, together 
with Vd equal to the identity matrix, gives a CKM matrix with mixings between the first 
and second generations and vanishing mixings between the rest. From a pragmatic 
point of view, it is important to remember that the CKM matrix is sensitive to kinetic 
terms and our computations are based on the leading order Kahler potential. In general, 
the Kahler potential receives corrections from various sources and computing extremely 
small quantities such as the second-third and first-third generation mixings without 
incorporating the corrections to the Kahler potential might be unreliable. Thus, while 
working with the leading order Kahler potential a reasonable goal might be to find models 
where the only non-vanishing mixing is between first and second generations. 



5.2 Up and down quarks both as D3D3 states 

Models constructed in this class are GUT models with multiple gauge groups such as 
left-right symmetric, Pati-Salam and trinification 0[8]. In these models the breaking of 

9 A realistic value of the second-third generation mixing can arise if some fields take vevs larger than 
the cutoff. However, the presence of such vevs most likely leads to the break down of the effective field 
theory. 
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the gauge symmetries occurs in two steps: first down to the Standard Model and then 
electroweak symmetry breaking at a lower scale. In the case of supersymmetric models it 
is typical that the breaking to the Standard Model preserves supersymmetry. As a result, 
some fields acquire supersymmetric masses at a very high scale and are non-dynamical 
in the low energy theory. 

Since both up and down quarks arise from within the singularity, the form of their 
Yukawa matrices is similar, although characterised by vevs of different fields. In the 
process of breaking the GUT symmetry, some fields can acquire a vanishing vev thereby 
giving the Yukawa matrices for the up and down quarks different structures. In order to 
illustrate this we begin our discussion with a left-right symmetric model constructed on 
the dP\ quiver. 

5.2.1 Left-right model 

A simple extension of the Standard Model is the left-right model with gauge symmetry 
SU(3) C x SU(2)l x SU(2)r x U(1)b-l- Its embedding in local models has been studied 
in [71181162] . We consider the left-right model from realised on dP\ |8], shown in Figure I3T1 
Prior to breaking SU(2)r the Yukawa matrices for the up and down quarks are both of 
the form of the dP\ Yukawa matrix, as in Eq. HI Each of the Yukawa matrices involves 
different Higgs fields. 




Figure 31: The dP\ Left- Right model. We vev the fields $ 6 i, A, A and B. 
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We initially break £7(2)^ x U(2) 2 R to £7(1) x £7(1) (generated by Tj + Q l ) by veving the 
field $6i. 

' 
v 
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The Higgs fields X 12 , Z 12 , and Y$2, are decomposed into 
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(34) 



where the last two entries in the row vectors denote the £7(1) quantum numbers. The 
right handed quarks transforming as (2r, 3) are broken into up and down quarks 
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where the first two charges indicate again the £7(1) charges. The superpotential for the 
quarks is 
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where if is the remaining field of $6i transforming as (2, —2) under the unbroken £7(1) 
factors. As discussed earlier the Yukawa matrix for the up and down quarks have the 
same structure, although the Higgs fields appearing in each matrix are different. 

We must ensure that both the abelian and non-abelian D-term equations are satified. 
The former can always be satisfied by tuning the FI terms, corresponding to resolving the 
singularity. The non-abelian D-terms require more fields to acquire vevs. A consistent 
solution can be obtained by giving the fields A, A and B vevs. These fields are shown 
in red in Figure [3TJ In Appendix [B] we show explicitly how the non-abelian D-term 
equations are satisfied. 

The effect of vevving D7D3 states is that we induce large supersymmetric mass terms 
which then are integrated out, rendering certain couplings absent at low energies. Any 
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terms involving only the high scale vevs may safely be integrated out. However, terms of 
magnitude v/A should be kept in the low energy superpotential. With the vev structure 
from the Appendix [Bj we find the following term in the superpotential which includes v: 



W D vX{ 2 B . 



This generates a large mass for Xf 2 , and the F-term for B requires that Xf 2 
at low energies the superpotential for the down quarks is 
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0. Thus 

(39) 



One can also find consistent vevs for the other Higgs fields which set the fields Zf 2 and 
Yq 2 to zero. The up quark Yukawa matrix is unchanged. 



5.2.2 Flavour mixing 

For the down quark Yukawa matrix in Eq. 
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the matrix Vh is 
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(40) 



where a 



^2*61 



and b 2 = 1 + 



AY* 



. Of course, this matrix could easily 



be obtained from the matrix in Eq. [26] by simply setting Xf 2 = 0. Similarly, by setting 
the fields Zf 2 and Yf 2 to zero one obtains the unitary matrices 
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where c 2 = 1 + 



AY* 



, d 
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Now let us compute the CKM matrix for Vd in Eq. HOJ Keeping in mind that the 
mixing angles in the CKM matrix are small, we take the limit 



Zf 2 $61 



(42) 



We will use the scaling of the CKM matrix in Eq. [2H to set the value of m later on. 

To obtain a CKM matrix which is close to the identity matrix, at zeroth order in e, 
V u must be the inverse of \Q. Motivated by this we take the ratios of fields in the Y u 
matrix to scale as 



vu 
A 12 

yu 
^12 



01 



r^ e nv* and 

Zf 2 A 



1JJ, 



(43) 



42 



Then to leading order in e 
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The resulting CKM matrix is 
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(45) 



One can see that it is not possible to obtain a realistic second-third generation coupling 
as the above form of the matrix fixes m = 1, n xz = 3 and therefore one of the second- 
third generation couplings is 0(e A ). One also can check that the other structures for Vd 
in (l4"T|) do not give a realistic CKM matrix. 



Moving to dP2 

This problem can be resolved by considering the model at the <iP 2 singularity, as de- 
scribed in Appendix After integrating the Higgs field Xf 2 and taking the limit 



e m < 1 



(46) 



the matrix Vh takes the form 
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Given this, for the CKM matrix to be close to the identity, V u at zeroth order in e must 
be of the form 
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This form is obtained by taking the scalings 
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As introduced earlier, we parametrise the scalings of these ratios by 
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All contributions which include ^42/ A are subleading. The matrix V u then scales as 
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With the following choice of scaling 

m = l,n x = l,n z = l,nq> = 2 
we obtain the correct hierarchical structure for the CKM matrix 
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6 CP Violation 

We have so far neglected the complex phase of the CKM matrix which is the source of 
CP violation in the Standard Model. A unique measure of the amount of CP violation 
was introduced by Jarlskog [44,45j. She introduced the matrix 



iC=[M u ,M d } = [Y u YlY d Y}\, 



(56) 
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and showed that the measure for CP violation is given by J, defined as follows 

detC = -2TB J , (57) 

where T = (m 2 — m^)(mj — m 2 )(m 2 — m 2 ) and B = (ml — m 2 )(m 2 — m d )(m 2 — m 2 d ). In 
terms of components of the CKM matrix, J is given as 

J = ImiVnV^V^) . (58) 

In the parametrisation of Equation [25] one finds (16] 

J = C1C2C3S1S2S3 sin 5 . (59) 

We note that to obtain the scaling in Equation [241 requires 

si ~ e, s 2 ~ e 2 and s 3 ~ e 3 . 

This implies J f=s e 6 sin 5. Thus, any model which reproduces the magnitudes of flavour 
mixing in the CKM matrix automatically provides an e 6 suppression in the Jarlskog 
invariant. Experimentally J is measured as 



<-oio>< KT 5 ^e 6 - 5 . (60) 

Thus in such models an explanation of the strength of CP violation reduces to showing 
that sin 5 ~ 0.5. In the following we give expressions for J and sin 6 in terms of Higgs 
vevs. This provides a further constraint on our models, but we cannot determine it here. 

Models with D3D7 and D3D3 states 

For non-zero CP violation Eq. ISH requires that the down quark Yukawa matrix, arising 
from D3D7 states, cannot be proportional to the identity matrix; a rotation between 
generations is needed. Recall that in the dP\ model presented in Section 15.11 such a 
rotation between the strange and bottom quarks was necessary to reproduce the CKM 
matrix. We take the down quark Yukawa matrix to be the same as in Equation [271 
denoting the phase of Y 62 Z 12 &qi/ A as 5 find 

J = - ■ TT7 — fn m — nr ■ rsr - ^ sino . (oil 

2 |X 12 |2+|Z 12 |2 | r62 |2 + ^ ( | Xl2 |2 + | Zl2 |2) ^ > 

Using the scalings of Section 5.1 one can indeed verify that the coefficient of sin 5 scales 
as e 6 

J^e 2 -e 4 sin5 = e 6 sin5 (62) 

in keeping with the general arguments presented earlier in this section. The correct 
magnitude of the Jarlskog invariant is obtained when the phase of Y62-^i2$6i/A is ap- 
proximately 0.5. 
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Models purely from D3D3 states 

In this section we discuss CP violation for the left-right symmetric model constructed in 
Section 5.2. On (IP2 with X d integrated out 

j _ \^ d \ 2 \X U \ 2 Im(Y u Y d Z d Z u $ u $ d ) ^ 
A 4 (m 2 . — m 2 d )(m 2 — m 2 a )(m 2 — m 2 ) 

where we have used the scalings for the ratios of fields in equations 06] and HH1 and 5 is 
the phase of Y u Y d Z d Z u $ u $ d . 

For both the D3D3 and D3D7 models we find that scalings introduced to achieve 
the correct hierarchies in the CKM matrix are consistent with the observed levels of CP 
violation if sin 5 ~ e 5 . 



7 Conclusion and outlook 

We have explored the phenomenological aspects of a well defined class of string models 
based on D-branes at toric singularities and substantially extended previous analyses 
that mostly concentrated on orbifold and del Pezzo singularities to an infinite class of 
models. This subject has been widely studied in the past few years from more theoretical 
perspectives and we have used the recently developed techniques in this field to extract 
phenomenological information. We have found several encouraging results. The fact 
that the number of families is bounded by the experimentally known case is remarkable, 
especially because in other constructions this number can be arbitrarily large. Also 
having a natural hierarchy of quark masses with one vanishing mass eigenvalue is very 
encouraging. In toric singularities the mass hierarchy at tree level purely arises from 
the presence of higher dimensional couplings. Instanton effects via E3 branes [63j can 
only contribute in very special circumstances. This contrasts with the intersecting brane 
models [64H66] . where such effects were necessary to create the hierarchy. 

The CKM matrix is more model dependent. At tree- level we can obtain the correct 
result as a unit matrix plus small corrections. These corrections can be affected from yet 
undetermined contributions from kinetic terms and therefore we cannot be conclusive 
before exploring them. If they turn out to be subleading, our results can be used to 
actually select some of the singularities over the rest. In this regard, the dP\ singularity 
stands out by the fact that it allows mixings among all generations for models with 
quarks from D3D7 and D3D3 states. The dP 2 singularity is the simplest example with 
all mixings for models with quarks purely from D3D3 states. In the limit of correct flavour 
mixing we automatically find the expected range of CP violation. It is interesting that 
both the mass hierarchy and the flavour mixing are achieved at the same time and are 
explained by demanding only the size of some ratios of Higgs vevs. 
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This work opens up many new vistas for exploration. As we mentioned, many of the 
results depend on several higgsings. A proper determination of the vevs of these fields 
by extremising the corresponding scalar potential should be done. Further phenomeno- 
logical issues such as neutrino masses and proton stability are more model dependent 
and we plan to study them in detail for particular realistic models. 

An important aspect of these models is the fact that the superpotential W (and 
therefore the Yukawa couplings) does not depend on the complex structure moduli. It 
is well known that in general W cannot depend on Kahler moduli due to Peccei-Quinn 
symmetries. But there are no such symmetries for the complex structure moduli. Never- 
theless toric singularities represent a large class of models where the complex structure 
moduli do not couple directly to the flavour dependant Yukawa couplings. This is rel- 
evant for supersymmetry breaking. Even in the cases where complex structure moduli 
could contribute significantly to supersymmetry breaking |6j, the soft breaking terms 
can still be flavour universal. 

Clearly the incorporation of these local models into fully fledged Calabi-Yau com- 
pactifications with the right singularity structure and all tadpole cancellation constraints 
satisfied needs to be addressed in order to have a proper string model. Also, to date, 
the techniques for non-toric singularities are not developed to the same extent as for 
toric singularities. For example, there is no dimer interpretation of non-toric singular- 
ities. Extension of these techniques to non-toric cases, such as higher order del Pezzo 
singularities would be desirable. We hope to come back to some of these issues in future 
publications. 
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A More dimer operations 



Here we give some more details on dimer operations. 

• Removing additional crossings when creating 3 copies of (1,1) paths. 
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Figure 32: Creating 3 copies of (1, 1) paths. Left: Creating them individually generates ad- 
ditional crossings (blue circles highlight additional crossings). Right: Removing the additional 
crossings changes the structure of the dimer locally in the same way as creating a single (1,1) 
path. 



B D-term equations 



We choose the following VEVs: 



v 



in 



.4 



v 2 
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(64) 



The choice for B is so that only Xf 2 becomes massive in the superpotential. The non- 
abelian D-term equations are 

XijtjfrXli = Yijtj k Yfo , (65) 

where Xy denotes the incoming fields, and Y^ the outgoing ones at a given node of the 
quiver. The vevs chosen above satisfy these D-term equations. 

We do not solve the D-term equations for the D7 gauge groups, which to be satisfied 
require the inclusion of D7-D7 strings. These D7-D7 strings can also be used to generate 
mass terms for some unwanted D3D7 states. 



C Generalising models to higher toric del Pezzos 

All models on lower del Pezzo surfaces can be embedded into models on higher toric del 
Pezzo surfaces. This can be desirable in the context of flavour physics as discussed in 
Section 5. 
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Going from dP\ to dP 2 introduces one additional field, ^42, and one additional gauge 
group "4". To keep the same matter content as in dPi the additional gauge group "4" 
must be at least as large as gauge group "2". One can vev ^42 to preserve the gauge 
group which is a diagonal combination of "4" and "2". 

For example, the left-right model on dP% discussed in Section 15.21 can be generalised 
to dP 2 . A vev 

breaks the two U (2) l to the diagonal one. The quivers of both theories are depicted in 
Figure 1331 The same procedure can be applied to going from dP 2 to dP%. 




Figure 33: A generalisation of a left- right model from dP\ to dP-2- The additional field *4 2 is indicated 
in red. 



D Determining all cycles in the quiver diagram 

In this appendix we present a systematic method to count the number of gauge invariant 
operators of a given length in a quiver |J7]. Gauge invariant operators are given by closed 
cycles in the quiver. For a quiver with n nodes, the adjacency matrix is an n x n matrix 
whose entries are given by the number of arrows going from node % to node j. In the 
case of dPi the adjacency matrix is 
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0/ 



(67) 



Raising this matrix to the I th power, the diagonal entries (a )u give the number of gauge 
invariant operators of length I. Note that all gauge invariant operators of length greater 
than n can be obtained as a product of gauge invariant operators of length less than or 
equal to n. 
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